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Abstract 

The eigenvalue spectra of staggered fermions with an Adams and/or Hoelbling mass term 
are studied. The chiralities of the eigenmodes reflect whether the chirality linked to the 
unflavored approximate (75 x 1) or the flavored exact (75 x £5) staggered symmetry is 
considered, and which one of the RR, LR, RL, LL eigenmode definitions is used. In 
either case a sensitivity to the topological charge of the gauge background is found. We 
discuss how to remove the leading cut-off effects of these actions by means of a properly 
tuned improvement term and/or the overlap procedure. The combination of Symanzik 
improvement and link smearing radically improves the properties of the physical branch. 



1 Introduction 

The so-called fermion doubling problem in lattice field theory has been addressed in many ways; 
the two most popular approaches are known as Wilson fermions [l] and staggered fermions [2] , 
respectively. A problem being cured in several ways usually indicates that none of the solutions 
is completely satisfactory in all respects, and the lattice is no exception to this rule. 

Wilson fermions represent the spinor components as internal degrees of freedom (i.e. with 
explicit 7-matrices), such that on a lattice with A" grid points the Wilson formulation entails a 
matrix of size 4N C N x 4N C N, where N c is the number of colors. Susskind "staggered" fermions 
are based on the observation that the eigenvalue spectra of naive fermions on interacting back- 
grounds are 4-fold degenerate (i.e. not 16-fold in 4 space-time dimensions as one might have 
naively guessed), and hence a reduction to 1-component spinors is possible (at the price of dis- 
tributing the spinor degrees of freedom over space-time), such that one ends up with a staggered 
fermion matrix of size N C N x N C N that corresponds to 4 species in the continuum. 

Wilson fermions are convenient for their conceptual simplicity (there is a 1-to-l correspon- 
dence between lattice and continuum flavor). Their main technical disadvantage is the breaking 
of chiral symmetry through the dimension 5 operator ip(—A)ip that is added to lift 15 out of 
the 16 naive species to a mass of order 1/a, where A is the gauge covariant Laplacian and 
a is the lattice spacing. Staggered fermions offer the advantage of a truncated form of chiral 
symmetry and of more speedy simulations (the size of the matrix is smaller). Their main tech- 
nical disadvantage is that the taste symmetry (among the 4 species) is not exact and that the 
reduction to a single species proceeds in different ways for sea-quarks (those which come from 
the functional determinant) and valence-quarks (those which stem from interpolating fields). 
This renders the rooted staggered formulation (at finite a) non-local and/or non-unitary which, 



in principle, could affect the universality class of the theory (though there is plenty of analytical 
and numerical evidence that this does not happen), see e.g. the review (3]. 

In either approach the main technical disadvantage can be mitigated. Replacing the gauge 
links Up(x) of the covariant derivative (V Al 0)(x) = [U /J ,(x)<p(x+jj l ) — UUx— jX)<f)(x— p)]/2 in the 
operator by smeared links V^x) (see below for details) reduces the amount of chiral symmetry 
breaking with Wilson fermions (particularly effectively when combined with a clover term) 
and of taste symmetry breaking with staggered fermions (see [3] for a guide to the literature). 
In the Wilson case chiral symmetry breaking can be completely removed through the overlap 
procedure (4j|5); unfortunately this increases the computational requirements by a factor 0(100). 

Recently, two modifications of the staggered action S$ = x[-Ds + m]x with D s = r/^V^ 
and rj^x) given below were proposed which go by the somewhat confusing name "staggered 
Wilson fermions". In this article we call them "taste-split staggered actions". In essence 
the proposal is to replace/augment the usual staggered mass term m(lcg)l) in Sg, where the 
notation is spinor®taste (6H8J , by taste non-singlet mass terms oc (1<8>£) which are designed 
to (partly or fully) lift the staggered near-degeneracy, i.e. some species get masses 2/ a or 4/a 
such that they decouple in the continuum limit. Such non-standard staggered mass terms 
were first considered in [8] where it was noticed that they break the remnant chiral symmetry 
Ds(x,y) — > e 16 ^ 9 D$(x, y)e le ^ 9 of the massless staggered action [with e(x) defined below] more 
severely than the usual mass term does. The Adams proposal [9] 

S A = xiv i y fl + r(M A + l)} X (1) 

builds on a mass term Ma — 1<S> £5, such that in the naive continuum limit two species stay 
massless, while two have mass 2r/a (where < r < 2 is akin to the Wilson parameter); the 
precise form of Ma will be given below. Similarly, Hoelbling proposes the two operators [To] 



S H ori = + r(M A + M Ho ri+2)]x (2) 

Sfisym = xln^f* + r(M Hs ym + 2)]x (3) 

which fully lift the staggered near-degeneracy, i.e. one species stays massless, while two receive 
a mass 2r/a and the last one is brought up to 4r/a. The two versions differ on how they break 
the rotational symmetry group R^ u of the massless staggered action; the precise form of M Ho ri 
and MH S y m will be given below. It is straight forward to write down the linear combination 

r 1 r i 

Sfinew = XiVt^fM + -(M A + M H sym + 3)]x = - [S A + SfisymJ (4) 

which lifts all three non-continuum species to the same doubler point 2r/a. 

Evidently, the attractive feature of these operators is that the matrices are just of size 
N C N x N C N, and still only 1 or 2 species survive in the continuum. The question has been raised 
whether the remaining symmetries of ([IJ or ([2j |3j [4]) are sufficient for taking the continuum 
limit (at a fixed pion mass) by tuning only the standard (relevant) mass term or whether 



other (relevant, marginal or irrelevant) operators need to be tuned simultaneously 11 17 . In 
this paper we investigate the eigenvalue spectra of these operators on interacting backgrounds 
and the chiralities (with respect to 75 (8)1, l®£s and 75®£s) of their eigenmodes. The idea 
behind is that the willingness or reluctance of the "bellies" to open up (and hence separate the 
branches) is indicative of how severe the fermionic operator mixing is - in the Wilson case both 
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link smearing and including the clover term tend to clear the first eigenvalue belly, since they 
suppress mixing between the dimension 5 Laplacian and the dimension 3 mass operator. 

In the proposals [9j[l0] (and in 11 16)) a notation is used which obscures the link to the 



original work j8]. This is why Sec.|2| contains a review of taste non-singlet mass terms using 
standard staggered notation. The core of the investigation, a look at the staggered eigenvalues 
and chiralities, is presented in Sec.[3j A brief discussion of some of the peculiarities of the 
Symanzik improvement program, when applied to these partly or fully undoubled staggered 
actions, is given in Sec.|4j Results where these taste-split staggered actions serve as kernel to 
the overlap procedure are shown in Sec. [5] A comparison to the eigenvalue spectra of Wilson- 
type actions is found in Sec.[6j Finally, Sec. [7] contains a summary. 

2 Review of staggered mass terms 

Define the Golterman-Smit staggered phase factors [8] 

Vlt (x) = (-1)£,<m^ , Ux) = (-1)£m<^ (5) 

which multiply to give 

V5 ( x ) = VlV2mV4 = (-l)«+*3 , Cs ( x ) = ClCaC 3 C4 = (_!)«»+*« ( 6 ) 

respectively. The product of the latter two phase factors reads (throughout the spinor®taste 
identification holds up to 0(a) corrections 



e(x) = rfeCs = (-l)^+*3+*4 < > 75 ® £5 (7) 

and induces the C/(l) e -symmetry and e-hermiticity e(x)D s (x,y)e(y) = D s (x,yY of the massless 
staggered action. After being lifted to a site-diagonal operator 

e(x, y) = T 55 (x, y) = {-1)^ +M ^8 x>y (8) 

it couples to the 75® £5 "Goldstone" state when used as an interpolating field [8). 
The (7 M <8>1) and (75(8)1) "taste singlet" operators are defined by 

rVo(z,2/) = r^(z,y) = ^r]^x)[U fl (x)5 x+flty + Ul(x-fi)5 x - fl)y \ (9) 

r 5 o(x,y) = r 6 (ar,j/) = — epermrir 2 r 3 r4 (10) 

perm 

and the (l®^) and (l®£s) "spinor singlet" operators are defined by 

Foti(x,y) = £n(x,y) = -(^{x)\pv(x)8 x +fr ty + Ul(x-jl)8 x -^y 
r 05 (x,y) = E 5 (x,y) = ■— ^ Sp^E^E^ (12) 

perm 

with the consequence that both r 50 and T 5 are 4-hop operators. Furthermore, the latter two 
operators relate to each other by a simple r 55 operation (either from the left or from the right). 
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In practice it can be advantageous to introduce the commutators in spinor and taste space 



r M ,(x,y) = ^(rvrv-rjv) <— ► 7^®i (13) 

Znu(x,y) ee ^("^ ~ S " S m) < — ► ( 14 ) 
respectively, with 7 Mi , ee ^[7^,7^] known as and ee ^[^,^]. Based on these one finds 

r 5Q (x,y) ~ -^(ri 2 r 34 - ri 3 r 24 + ri 4 r 23 + r 23 ri4 - r 24 ri 3 + r 34 ri 2 ) (15) 

Fo5{ x ,U) — —q ^12^34 — ^13^24 + 23 + ^23^14 — "24^13 + ^34^12^ (16) 

where the near-equality is exact if no SU(3) backprojection in intermediate steps is applied. 
The main advantage of the form ( 15fl6 ) is that one can apply SU(3) backprojection twice, in 



each case the argument being a sum over products of just two color matrices. 

Following Golterman and Smit (8), everything which is a singlet in spinor space is called a 
"mass term" . This leads to the categorization of potential mass terms as 



1®1 


< — > 1 


(0-hop, taste-scalar) 


(17) 




<— ► ^ 


(1-hop, taste-vector) 


(18) 


1®£^ 


< — > 


(2-hop, taste-tensor) 


(19) 




< > 1= M =5 


(3-hop, taste-pseudovector) 


(20) 


1®£ 5 


< > =5 


(4-hop, taste-pseudoscalar) 


(21) 



in terms of hermitian operators. As was pointed out in [11], an acceptable mass term is 1<8>1, 
l®£/w) 1®^5 or a combination thereof. In this terminology the Adams term in (JIJ) is the 
taste-pseudoscalar mass r s = S 5 . And the Hoelbling terms in (J2j [3]) represent unsymmetrized 
and symmetrized versions of the taste-tensor mass E^, respectively, i.e. Mr-oh = S12 [to be 
used together with M A , see Q] and M Hs ym = (H12 + S 34 + S i3 - H 24 + H 14 + H 23 )/v / 3 [to be 
used alone, see ([3p]. The operator Q, finally, involves all three types of valid mass terms. 

Note that there is some freedom in the definition of the operators and of the taste 

non-singlet mass terms. For instance with the Adams operator as defined in Q and the choice 
Ma = 75® £5 it follows that the physical branch has £5 = —1. Hence, with this convention the 
actions of 75(8)! and 75(8)^5, in the physical branch, on a topological mode will differ in sign. 



3 Eigenvalues and chiralities 

The eigenvalue spectrum on a thermalized gauge background is a key to get an impression 
of how suitable a given fermion discretization is in practical terms. An issue with staggered 



fermions is their sensitivity to the topological charge of the gauge configuration 18-21 . 

The e-hermiticity ensures that the eigenvalues of the massless operator occur in pairs ±iA 
on the imaginary axis. There are no exact zero- modes of D$; instead 4\q\ would-be zero-modes 
show up. Both their separation from the non-topological modes and their chiralities (see below) 
depend on how much the staggered action is "improved" through link smearing. 
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Figure 1: Eigenvalue spectra of the four unimproved operators (csw = 0) at m = for up to 
three levels of HEX smearing - the smearing seems essential for the bellies to form. 



Fig. [I] displays the eigenvalue spectra of the operators (JT]^4|) on a 6 4 lattice (generated with 
= 5.8 Wilson glue) of unit topological charge, \q\ = 1. The fermion boundary conditions are 
periodic in all directions; with antiperiodic boundary conditions in the 4-direction mild changes 
would occur in regions of large (positive or negative) imaginary part. Throughout, we use gauge 



links V^(x) which have undergone 0, 1, or 3 levels of HEX smearing 22 . Backprojection of the 



mean of the n-hop paths to SU(3) could be applied, but in the present work this is not done. 

Our first observation in these plots is that without link-smearing the "bellies" have a hard 
time opening up (as was already observed in fll]), but the situation improves considerably 
upon applying 1 or 3 steps of smearing. All operators have a bare (taste-scalar) mass m — 0, 
but their physical (left-most) branches cross the real axis at distinctively non-zero values. This 
is a sign of additive mass renormalization, and the plots show that the smearing reduces this 
effect significantly. Last but not least the physical branches have one [or two for 0, since 
this operator is doubled] exactly real modes which tend to get "soaked into the belly" even for 
the higher smearing levels. And of course all spectra are symmetric with respect to the real 
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4D: NC=3, 0=5.8, L=6, T=6, |q|=1 , c_SW=0 4D: NC=3, p=5.8, L=6, T=6, |q|=1 , c„SW=0 




Figure 2: Needle plots of the four L/R-chiralities of the unimproved Adams operator (csw = 0) 
with respect to r 50 (green circles) or r 55 = e (blue stars), with three HEX smearings. 



axis, which implies that the determinant is real (this is a consequence of the requirement of 
e-hermiticity which ruled out two possible mass terms, cf. Sec.[2|. 

The next step is to look at the chiralities of the pertinent eigenmodes, i.e. the values of the 
operators r 50 , r 5 and r 55 when sandwiched between the staggered eigenmodes. The standard 
operator D$ or D$ + m is a normal operator, i.e. [D$, D$] = 0, and only r 50 is sensitive to the 
topological charge of the background (see below). The taste-split staggered operators (fil-ffl are 
non-normal operators, i.e. [D, D*] 7^ 0. This implies that the left-eigenvectors, which satisfy 
(Li\D = (Li\Xi, are not just the hermitian conjugates of the right-eigenvectors, which satisfy 
D\Ri) = \i\Ri), though they share the eigenvalue Aj. The situation is now analogous to that 
of the Wilson operator which is also non- normal [23]. This implies that for a given sandwich 
operator T there are four chiralities, (L;|r|Lj), (Li\T\Ri), (Ri\T\Li), (Ri\T\Ri), where |L») is 
the hermitian conjugate of the (normalized) left-eigenvector (L»|, and (Ri\ is the hermitian 
conjugate of the (normalized) right-eigenvector \Ri), with % running over all modes. 

Fig.|2]displays these four options with the sandwich operators r 50 and r 55 for the eigenmodes 
of the Adams operator uty. The chiralities are plotted in the ^-direction above the respective 
eigenvalue A = x + iy in the complex plane. The four options collapse effectively into two, 
as it must be, due to the e-hermiticity of the operator. The (L\T\L) or (R\T\R) option shows 
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4D: NC=3, p=5.8, L=6, T=6, |q|=1 



4D: NC=3, p=5.8, L=6, T=6, |q|=1, C_SW=0 
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Figure 3: Chiralities of the standard staggered action and of the unimproved Wilson operator. 
For better visibility the results for (.|r 50 |.) and (. |r s I -) are displaced in the real direction. 



4D: NC=3, p=5.8, L=6, T=6, |q|=1, C_SW=0 
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4D: NC=3, p=5.8, L=6, T=6, |q|=1, C_SW=0 



4D: NC=3, p=5.8, L=6, T=6, |q|=1, C_SW=0 
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Figure 4: Needle plots of the chiralities with respect to r 55 of the unimproved operators (csw : 
0) above their eigenvalues. The LR-dehnition of (.|r 55 |.) and three HEX smearings are used. 



two distinctively non-zero physical modes (on a configuration with \q\ = 1), depicted at the 
position of the two exactly real eigenvalues, both for T = T 50 and T = r 55 (peeking into 



7 



the +z and — z directions, respectively). With the (L|r|i?) or (R\T\L) option the chiralities 
tend to be even more pronounced, and in either case it holds that the magnitude of (.|r 55 |.) 
exceeds the magnitude of (.|r 50 |.). Finally there is an equal number of chiral modes in the 
unphysical branch (near x — 2) with identical orientation for r 50 but opposite orientation for 
r 55 (in accordance with Karsten-Smit [18]). The most surprising lesson is that both r 50 and 
r 55 prove sensitive to the topology of the gauge background; this happens in sharp contrast to 
the standard staggered case (cf. next paragraph). And the recommendation of 23 that for any 
T one should focus on the {L\T\R) option continues to be valid. 

Fig. [3] displays - for comparison - the chiralities of the standard staggered action and of the 
Wilson action on the same \q\ = 1 gauge background (only half of the spectrum is shown in 
the latter case). In the staggered case there is only one type of eigenmode, due to (Vi\ = \Vi)^ , 
while in the Wilson case only (Li\^\Ri) is shown. In the staggered case (Vi|r 55 |Vi) is absolutely 
flat (as required by {e,D s } = 0), while (^|r 50 |V^) (plotted at x — 1) has 4 modes which peek 
upwards, and (Vi|r 5 1 V^) (plotted at x — 2) is wiggly but not sensitive to topology. The Wilson 
operator shows the expected behavior - one chiral mode in the physical branch, four oppositely 
oriented modes near x = 2, six originally oriented modes near x = 4, four oppositely oriented 
modes near x = 6 (not shown) and one originally oriented mode near x = 8 (not shown). 

Fig. [J] displays the (L|r 55 |i?) chiralities for the four operators ([IfQ above their eigenvalues. 
The Adams operator has two chiral modes in the physical branch and two unphysical modes 
with opposite orientation. The unsymmetrized Hoelbling operator has one chiral mode which is 
physical, two modes of unequal orientation near x = 2 and one oppositely oriented mode near 
x = 4. The symmetrized version has one chiral mode which is physical, two oppositely oriented 
modes near x = 2 and one mode with original orientation near x = 4. The new operator Q, 
finally, has one chiral mode which is physical and three unphysical modes near x = 2 (two of 
which are antiparallel, one of which is parallel to the chirality of the physical mode). In each 
case it was checked that the choice r 50 (instead of r 55 ) or the naive option (R\.\R) (instead of 
(L\.\R)) bring exactly the kind of change that one would anticipate from Fig.[2j 



4 Results with Symanzik improvement 

So far the investigation of the eigenvalues and eigenmode chiralities of the taste-split staggered 
actions ([l}j4]) leaves us with the impression that they are technically rather close to the usual 
Wilson action - with chiral symmetry breaking and all the consequences of non-normality, 
e.g. (L\ and \R) eigenmodes. This raises the question whether some of the remedies which 
have proven useful in taming these effects with the Wilson action might be taken over and/or 
adapted to these taste-split staggered actions. The successful remedies were link smearing, 
clover improvement and the overlap procedure - they all mitigate (in the last case eliminate) 
the effects of chiral symmetry breaking, and they can be used in various combinations. 

That link smearing proves useful for staggered actions with taste non-singlet mass terms 
has already been shown in the previous section - in Fig.[T]a clear improvement of the overall 
properties of the eigenvalue spectra with an increased number of smearing steps is evident (and 
consequently in Figs.[2j |3j [5] only the results for three smearing steps were displayed). 

The next step is Symanzik improvement [24] . In principle this is a program to be carried out 
in strict analogy to the Wilson case - one writes down a complete list of operators, ordered by 
their mass dimensions, and eliminates those which are forbidden by symmetries or redundant 
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Figure 5: Eigenvalue spectra of the four improved operators (csw = l) at m = for up to three 
levels of HEX smearing - the smearing seems essential for the bellies to form. 



by the equation of motion. In the Wilson case one finds that only the d = 5 term ip^^F^i/j 
needs to be added to the action [25], with a tuned coefficient csw/2 in front to cancel all O(a) 
cut-off effects (thus eliminating all signs of chiral symmetry breaking to this order) [26|[27| . 

With staggered fermions the situation is a bit more delicate, due to the "staggered" rep- 
resentation of the spinor degrees of freedom. The point is that the operators = r/^C^ and 
S M = C/iC/i formally leave the mass dimension d of a fermion bilinear invariant, but they do 
involve a fermion hop [made gauge invariant through a link U^x) or V^(x)]. In other words, the 
strict coincidence of the mass dimension of the operator (as is relevant for the Symanzik anal- 
ysis) and the maximum number of hops (that affects the strength of the mixing or the "noise" 
of the operator) is now broken. It seems that this requires a more in-depth analysis |l7j. Here 
we give an incomplete and redundant list of staggered bilinears (with "+h.c." implicit) 

d = 3-. ^{i,r M ,s Ml i^,r M s M ,sJ,...}v (22) 

d = 4 : H^D^ r MJ D M , S„£>„, (23) 
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SW=1 




Figure 6: Needle plots of the LL (left) and LR (right) chiralities of the Adams operator |ip 
(top) and of the new operator (H) (bottom), with respect to T 50 (green circles) or r 55 = e (blue 
stars), with Symanzik improvement (cgw = l) and three HEX smearings. 



d = 5: ij{D% V ,Dl, C,D% T,D% ~,D% T, U F, U , E, U F, U , ...}^ (24) 

along with a statement that most of them can be eliminated by means of the arguments 
mentioned above and the identification of operators which differ only by a factor (1®^). 
The reasoning sketched in App.|A| suggests that only a single operator needs to be included 

- ^{r 12 F 12 + ... + r 34 F 3 4 + F l2 r l2 + ... + F 34 r 34 } = -^{ £ rw + h.c.} (25) 



fl<V 



to be compared to — X^<i/ l\i.vF^x,y with Wilson fermions (in either case factors are chosen 
such that csw = 1 at tree level). Note that F^ v is site-diagonal, but not color-diagonal, whereas 



r MJy is neither site- nor color-diagonal. The symmetrization in (25) ensures gauge-covariance 
(in contradistinction to gauge-invariance in the Wilson case) and hermiticity. 

Fig. [5] displays the eigenvalues of the four operators ([l}4]) with the improvement term (25) 



on the same gauge configuration (with |g| = 1) as in Sec!|3 Relative to Fig.[T]one finds a clear 
amelioration of the behavior of the physical (leftmost) branch of eigenvalues; it is now much 
thinner and the section close to the origin is much more akin to a shifted Ginsparg- Wilson circle. 
In particular the tendency of the exactly real modes to get "soaked into" the belly is almost 
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4D: NC=3, p=5.8, L=6, T=6, |q|=1, C_SW=1 
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Figure 7: Needle plots of the LL (left) and LR (right) overlaps (.|r 00 |.) (yellow circles) and 
taste chiralities (.|r 5|.) (red stars) for the Adams operator |Ip (top) and for the new operator 
(pi) (bottom), with Symanzik improvement (csw = l) an d three HEX smearings. 



removed, except for the unsymmetrized operator (J5]). Perhaps the most remarkable feature is 
that for the unsmeared operators the improvement is hardly useful. It takes the combination 
of smearing and improvement to get a profound effect - in strict analogy to what is observed 



for the Wilson case 22 . For clarity we add that the links in the operator, in Y^ and in F^ v 
have all undergone the same kind of smearing. Finally, let us add that it was explicitly checked 
that the alternative operator J2fi<u E^i^ + h.c. does not lead to a similar improvement of the 
eigenvalues in the physical branch. 

Fig. [6] displays the same kind of stem plots for the r 50 and r 55 chiralities of the eigenmodes 
of the operators ([!]) and Q that were presented previously without improvement. Relative 
to Fig. [2] the most significant change is that the difference between (L\.\L) and (L\.\R) in the 
physical branch is now less pronounced. This is a hint that the non-normality of the operator, 
when restricted to the physical subspace, is reduced by the improvement term. 

Fig.[7]displays a feature that was not discussed in the unimproved case. The overlaps (L|1|L) 
and (L|l|i?) are shown for the operators ([T| and Q - the former ones are trivially one, but 
the latter ones indicate how normal the operator is, and indeed this figure is close to 1 in the 
physical branch. In addition the taste chiralities are displayed; both (L|r 5|L) and (L|r 5|-R) 
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Figure 8: Needle plots of the chiralities with respect to r 55 of the improved operators (csw — 1) 
above their eigenvalues. The LR-dehnition of (.|r 55 |.) and three HEX smearings are used. 



tend to assign each branch a fixed taste chirality which is not sensitive to topology of the 
background (as mentioned earlier, with our choices the taste chirality in the physical branch 
is £5 = —1). The conclusion is that the (L\.\R) taste chiralities in the interacting case look 
similar to the free-field results presented in [Ti] . 

Fig. [8] displays the (L\T^\R) chiralities of all four operators ([lJQ with the clover term (25). 
These chiralities are very well pronounced (i.e. very close to ±1) and linked to the topology of 
the gauge background; a cut at ±0.5 is well suited to identify the chiral modes. 



5 Results with overlap projection 

The definition of the massless overlap Dirac operator takes the form |4||5] 

D ov = ^(l + X(XtX)- 1 / 2 ) = ^(l + (XXt)-V2 X ) (26) 

where X = aD^ c — p and the kernel may be any undoubled fermion action [as is the case 
with (|2]-[4])] or a doubled one where the tastes in the physical branch share one chirality [as is 
the case with (HI)]. Evidently, the overlap inherits the multiplicity of the kernel operator, e.g. 
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4D: NC=3, p=5.8, L=6, T=6, |q|=1 , C_SW=1 4D: NC=3, 0=5.8, L=6, T=6, |q|=1 , C_SW=1 




Figure 9: Needle plots of the chiralities with respect to T 55 of the overlap actions with improved 
kernels (csw — !)• All L/R-definitions of (.|r 55 |.) are equivalent; three HEX smearings are used. 



two- fold with Q as kernel. The canonical value of the shift parameter 0<p<2isp=l, 
a choice which we shall adopt in the following. Because the eigenvalue spectra are somewhat 
boring (the spectra lie on the shifted unit circle), we proceed directly to the chiralities. 

Fig.[9]displays the (.|r 55 |.) chiralities of the overlap action built from the kernels ([ifQ. These 
operators are normal, i.e. the left-eigenvectors are daggered copies of the right-eigenvectors, viz. 
(Li\ = \Ri)', and the distinction between the four L/R-versions of (.^55).} is now obsolete. The 
overlap procedure being a projection, together with the hermiticity of the kernel, ensures 
that the new (.|r 55 |.) chiralities are exactly ±1 (topological modes) or exactly (non-topological 
modes). As a side-effect in the conglomerate of unphysical branches (near x — 2) only one type 
of chirality survives [contrary to what a naive shift of the corresponding modes of the kernel 
operator in Fig.[8] would suggest]. Thus the needle at z = (x,y) = (2,0) is two-fold populated 
with the Adams kernel, but one-fold in all other cases. 

Fig. 10 compares the (.|r 55 |.) chiralities of the overlap operators with taste-split staggered 
kernel to the (.|r 50 |.) chiralities (still on the same \q\ = 1 configuration). While the former are 
exactly ±1 (for the topological modes) the latter are not, in spite of the overlap projection. As 
discussed before this is a consequence of the kernels being hermitian but not hermitian. 

As a practical point let us remark that the versions with unimproved kernels (csw — 0) look 
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4D: NC=3, B=5.8, L=6, T=6, |q|=1 , C_SW=1 4D: NC=3, 0=5.8, L=6, T=6, |q|=1 , C_SW=1 



E 




Figure 10: Needle plots of the chiralities of the overlap actions based on the Adams kernel |Ip 
(left) or the new kernel (pD (right), with respect to T 50 (green circles) or r 55 = e (blue stars). 
Either kernel uses Symanzik improvement (csw = l) an d three HEX smearings. 



4D: NC=3, B=5.8, L=6, T=6, |q|=1 , c_SW=1 4D: NC=3, 6=5.8, L=6, T=6, |q|=1 , C_SW=1 




Figure 11: Needle plot of the (.|7s|.) chiralities of the eigenmodes of the overlap action with 
clover improved (csw = 1) Wilson kernel (left) and (R\T 05 \R) taste-chiralities of the overlap 
action with clover improved (csw = lj Adams kernel (right). 



superficially identical to the versions with clover improved kernels (csw = 1) that were presented 
in the last two figures. Given that the combination of clover improvement and link smearing 
facilitated the separation between the physical branch and the unphysical branch(es), it is 
hardly surprising that the condition number in the overlap construction is significantly smaller 
in the latter case (which reduces the order of the polynomial/rational approximation to the 
sign function and thus the computational requirements). In short, it is strongly recommended 
to equip the kernel operator with a clover term and overall link smearing. 

Fig. 11 contains two addenda to the overlap theme. In the left panel the (.|7s|.) chiralities 
of the overlap operator with Wilson kernel are shown, and the analogy to Fig. 10 is evident. 
In the right panel the taste chiralities (.|r 5|.) of the overlap operator with Adams kernel are 
shown, and the smooth pattern that would show up in the kernel only in the (L|.|i2) version 
(cf. Fig. [7]) is now generic (as the L/R distinction is gone). 
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Figure 12: Eigenvalue spectra of the operators ([Ip and ([4]) and of the Wilson and Brillouin 
operators on a gauge configuration with topological charge \q\ = 1, with csw = (left) and 
csw = 1 (right), after 1 (top) and 3 (bottom) steps of HEX smearing. 



6 Eigenvalue comparison to Wilson-type actions 



Since the taste-split staggered operators (|lj-|4j) were found to behave, on many practical issues, 
like Wilson fermions, it seems instructive to compare their eigenvalue spectra (still on the same 
\q\ = 1 configuration) to those of the Wilson [I] and Brillouin [28] action. 



Fig. 12 gives such a comparison, without improvement (left) and with tree-level clover im- 
provement (right). Throughout gauge links V^x) are used which have undergone 1 or 3 steps of 
HEX smearing [22] . The striking feature with csw — is that the lowest non-topological modes 
of all operators sit essentially in the same place (with two nearly degenerate copies in case of the 
Adams operator), but the tendency of the topological modes to get "soaked into the belly" is 
different (the Adams operator fares best, the Wilson operator fares worst on this point). With 
csw — 1 the situation is just opposite; the "soak in" phenomenon seems cured with all opera- 
tors, but now the position of the physical branch varies, and this time the Brillouin operator 
performs best. A peculiar issue with the latter action is the "thorn" of eigenvalues that grows 
from the doubler point 2 into the belly, once the clover term starts spreading the 15 unphysical 
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modes. It is not clear whether this is an advantage (it might facilitate topology changes) or a 
disadvantage (topology might be less clearly defined) of this action - in the latter case it might 
be cured by multiplying its clover term with a factor (1 + A Bri /4), with symmetrization, where 



A 1 is defined in 128 1. In short these spectra underline the value of the improvement term (25) 



and suggest that - at least for some applications - the actions |4j) might fare reasonably well. 

7 Summary 

The main findings of this paper may be summarized as follows: 

1. The discretizations ([I]) and ([3]), as proposed by Adams [9] and Hoelbling 10 , respectively, 
and the new linear combination (2J yield staggered fermions with non-standard (i.e. taste 
non-singlet) mass terms (cf. Sec. 21). They are "hybrids" in the sense that they distribute 
the spinor degrees of freedom over spacetime in the manner of staggered fermions, while 
technically being close to Wilson fermions, with additive mass renormalization and the 
same consequences of non- normality ((L\ and \R) eigenmodes). 

2. The main technical challenge of staggered fermions - the non-coincidence of the exact 
(tasteful) e = r 55 and the non-exact (taste-singlet) r 50 chiral symmetries - persists. 
What changes is the details of how the operator is sensitive to the topological charge q 
of the gauge background. With just a scalar mass term, ( V"| r 55 1 V^) is zero between all 
eigenmodes of the operator and only (V^l^ol V) is sensitive to topology. With a sufficiently 
large taste-tensor or taste-pseudoscalar mass term, as is the case in (|T]^4|) with r = 1, the 
matrix element (L|r 50 |-R) continues to be sensitive to topology (with equal sign in the 
physical and the unphysical branches), while (L\T 55 \R) acquires an even better sensitivity 
to topology (with net zero sign across all branches). 



The physical branches of the operators (|lj-|4j) show a cancellation-free sensitivity to topol- 
ogy with 2\q\ or \q\ exactly real modes. But they also suffer from the same symptoms 
of chiral symmetry breaking with induced 0(a) cut-off effects as Wilson fermions, and 
it seems thus natural to try similar remedies. In this article it is conjectured that the 



leading Symanzik improvement term to these actions takes the form (25), and the con- 
nection to the overlap procedure suggests that the tree-level improvement coefficient is 
again csw — 1- The combination of Symanzik improvement and link-smearing is found 
to have a pronounced (and beneficial) effect on the formation of a "belly" between the 
physical and the unphysical branches in the eigenvalue plot (which, in turn, is seen as a 
sign that the mixing with higher-dimensional fermionic operators is suppressed). 

4. The kernel operators (jl]) and Q, once equipped with link smearing and a clover term, 
bear the promise of a cheap overlap construction (in contrast to the mild savings that 
were found without these ingredients [II]) due to a moderate spectral range over which 
the sign function is to be constructed. Accordingly, a combination of these two overlap 
operators might be an attractive option for simulating QCD with 2 + 1 + 1 active flavors. 



Acknowledgments: The author wishes to thank Christian Hoelbling and Stefan Krieg for 
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A Leading terms generated by the overlap procedure 



To facilitate the analysis of the necessary counterterms for 0(a) improvement of the taste-split 
staggered actions ([TfQ let us begin with a reflection on the situation with Wilson fermions. The 
idea that for Wilson fermions the overlap procedure of [ij|5] would automatically generate the 
necessary improvement term is mentioned in the review by Niedermayer p9j. Unfortunately, 
this account is not very verbose, and this is why we attempt a short summary. 
One essential ingredient in what follows is the continuum relation 

If) 2 = D^D vlv = (^{7m,7,} + \[l^l v ])D»D v = D 2 + J2 (27) 

with o^y = |[7 M , 7^] which is also dubbed 7 M „ in Sec.|2j The other ingredient is the definition 
(26) of the overlap procedure with X = aD^ e — p. This definition does not rely on any special 
property of the kernel operator D^ e . With an argument which is 75 or e hermitian the procedure 
can be recast in a form which involves the sign function; this will be relevant for numerical 
applications with as a kernel but not for the analytical considerations below. 

With the Wilson operator as kernel one starts from the transcription (here we follow [30]) 



Av = EKv,-^A,} ~ P- a -D 2 + 0(a 2 ) (2* 

A 1 



in terms of continuum operators. Here V M and A M denote the gauge covariant first and second 
discrete derivative, respectively, the latter one to be distinguished from the continuum D 2 . 
This implies X = -p + ap - °^D 2 + 0(a 3 ) and X* = -p - afl - ^D 2 + 0(a 3 ) and thus 

X*X = p 2 - (l-p)a 2 D 2 -a 2 Y, + 0(a 3 ) (29) 

in the physical branch. Next we use the expansion (1 + z)~ x l 2 — 1 — |z + 0(z 2 ) to arrive at 
P^)- 1/2 = ~ (l + ^a 2 D 2 + ^ £ a, v F, v + 0(a 3 )) (30) 



and upon multiplying this with X we find 



2 2 
a ^ a „ 9 a 



X(XtX)- 1 / 2 = -1 + -J) - ^- 2 D 2 - ^ £ a, v F, v + 0(a 3 ) (31) 



p^ 2p 2 2p 2 



with the consequence that the overlap operator relates to the continuum Dirac operator like 

D ov = p-^D 2 -^Y, <V*^ + °( fl2 ) • ( 32 ) 

^P ^P n<v 

The second term may be removed by a field rotation, but the third term indicates that 0(a) 
improvement of Dw calls for a term — Y^^c^F^ with csw = 1 at tree level. 

The task is now to pipe the taste-split staggered operators ([l}j4]) through the overlap pro- 
cedure and to see which counterterms are generated. First we need to transcribe the staggered 
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operator in terms of continuum operators on the blocked lattice (b = 2a). For ([I]) we find 



= E{(7^1)V M - 9 (75®^ 5 )A4+r(l®^ 5 ) + r(l®l) 



)75 D 2 + 0(6 2 ) 



(33) 



where in the second line [which is supposed to capture the effect on the physical branch in terms 
of continuum operators] we use that the violation of Lorentz symmetry takes place exclusively 



in the taste space. Hence, the only relevant difference to (28) is an additional factor of 75. This 



implies Y = -p + bp - y 75 £> 2 + 0(b 3 ) and yt = -p - bp - \^D 2 + 0(b 3 ) and thus 

yty = p 2 _ (i_ p75 ) 6 2 jD 2 _ b 2 ^ 7 ^ + 0(fe3) 

fl<v 



(34) 



in the physical branch. With the same expansion as before we obtain 
(yty)-V2 



P 



and upon multiplying this with Y we find 



1 ' 1 + i?» !fl, 4 Ev "'" + * 



(35) 



y(yty)-V2 = _ 1+ °^ 



p ■ 2p 2 2p 2 



E7^ + 0(& 3 



(36) 



which agrees with the expression (31) which we found in the Wilson case. The bottom line 



is that in order to cancel the 0(b) effects we need to add a term which acts on the physical 
fi<ul^Ffiu, with csw = 1 at tree level. To improve the actual action (|T1) 



branch like -5sp£ 



the obvious replacement is 7^ — > T^. Finally, to ensure gauge covariance and hermiticity this 



expression needs to be symmetrized in the manner of (25) 



None of the manipulations listed above did refer to the details of the taste lifting term ex r 
in (flf|4|); in fact an essential ingredient was that this term is effectively in the physical branch. 



It follows that the operator (25) removes the leading cut-off effects in all taste-split staggered 
actions, with csw = 1 at tree-level, regardless of the multiplicity of the physical branch. 
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